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We show here that due to non-adiabatic couplings in decaying systems applying the same time-dependent 
protocol in the forward and reverse direction to the same mixed initial state leads to different final pure states. 
In particular, in laser driven molecular systems applying a specifically chosen positively chirped laser pulse or 
an equivalent negatively chirped laser pulse yields entirely different final vibrational states. This phenomenon 
occurs when the laser frequency and intensity are slowly varied around an exceptional point (EP) in the laser 
intensity and frequency parameter space where the non-hermitian spectrum of the problem is degenerate. The 
protocol implies that a positively chirped laser pulse traces a counter-clockwise loop in time in the laser param- 
eters' space whereas a negatively chirped pulse follows the same loop in the clockwise direction. According 
to this protocol one can choose the final pure state from any initial state. The obtained results imply the intrin- 
sic non-adiabaticity of quantum transport around an EP, and offer a way to observe the EP experimentally in 
time-dependent quantum systems. 



For open quantum systems where the effective Hamiltonian 
is non-hermitian (NH) the non crossing rule (01 is replaced by 
an intersection of two complex energy levels associated with 
two eigenfunctions of the NH Hamiltonian that have the same 
symmetry. Let us consider 2x2 Hamiltonian matrix H which 
depends on potential parameters q\ and q2- These can be for 
instance the laser frequency and intensity when light interacts 
with two normal modes of a molecule. 

In open quantum system where the effective Hamiltonian 
is NH all matrix elements can attain complex values. The 
complex diagonal matrix elements are associated with meta- 
stable (resonance) states, such that —2ImHu and — 2imU22 
are the decay rates of the meta-stable states. The eigen- 
values of this NH Hamiltonian are degenerate when A = 
(Hn — H22) 2 + 4-ffi2i?2i = even though all matrix el- 
ements are different from zero. This situation is very differ- 
ent from the hermitian (standard) case where crossing requires 
fi\2 = H21 = and H11 = H22- At the crossing point a non- 
hermitian degeneracy (NHD) is obtained when the following 
two equations are satisfied: 



Re{Hn 
Im[H n 



H22] 
- H22] 



T2Im[H 12 H2i} 1/2 
= ±2Re[H 12 H2i} 1/2 . 
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(2) 



NHD is very different in its nature from hermitian degen- 
eracy. NHD is obtained at the crossing point denoted by 
(qf p , q2 P ), where the two eigenvalues coalesce and form a 
branch point (BP) in the complex energy spectrum. At the 
BP the first order derivatives of the eigenvalues with respect 
to qi or q2 acquire infinitely large values (see for example 
Chapter 9 in Ref.|2|). This BP is also known as an excep- 
tional point (EP) in the energy spectrum [Q 0]. Moreover, at 
the BP (EP) not only the eigenvalues coalesce but also the cor- 
responding eigenvectors. Such a phenomenon can never occur 
in standard QM. In NHQM as q x ->■ qf p and q 2 -> qf p the 
two bi-orthogonal eigenvectors of the complex non-hermitian 
Hamiltonian matrix coalesce. Rather than two different bi- 
orthogonal eigenvectors we get only one eigenvector which 



is self-orthogonal (with respect to the c-product)[|5l @]. As 
proved in Refs.flll, the NHD (i.e., EP) is a typical phe- 
nomenon in NHQM. 

Since there is no analog to this situation in hermitian QM 
it was believed for many years that the EP is a mathematical 
object only as it appeared mostly in complex scaled Hamilto- 
nians 0. Yet, even as a mathematical object it was found to 
be a useful concept that helps to explain experimental results 
which could not be explained otherwise (see Ref.[9J] where 
the mysterious sharp peaks in the cross section measure- 
ments of electron scattering from hydrogen molecule were ex- 
plained for the first time). In the last decade it became clear 
that EPs are not only mathematical objects but play a major 
role also in actual measurable phenomena. Different man- 
ifestations of the EPs have been described in optic s lfioll . in 
superconductors Ull, in q uantum phase transitions in a system 
of interacting bosons lll2ll . in electric field oscillations in mi- 
crowave cavities ll 1311 . and in PT-symmetric waveguides 114 1. 
So far there are no experimental results regarding EPs in 
atomic, molecular, or biophysical systems. 

Before proceeding we should mention the most striking 
phenomenon induced by EPs which has no equivalent in her- 
mitian QM: the state exchange phenomenon. The state ex- 
change phenomenon can be illustrated as follows: Consider 
an arbitrarily variation of the two parameters (91,(72) which 
depend on angle variable ip in a closed loop around an EP 
(qf p , qf p )- The two instantaneous eigenvalues are given by, 



E ± (<p) = 



ffu + H22 ± v'AM 



(3) 



where the quantity A(ip) makes a circle around the origin 
in complex plane with a change of <p. Thus, it is easy to 
see that E±(0) = E^(2n). Instead of the Berry phase 
which is obtained when cycling around a conical intersec- 
tion, when cycling around an EP one state flips into the other 
(see Chapter 9 in Ref.yfl). To the best of our knowledge 
the only measurement of this state exchange phenomenon 
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was carried out by Richter and his co-workers in microwave 
experiments 1 1311 . The association of the state exchange phe- 
nomenon with molecular system was made by Lefebvre and 
his co-workers |15 1 . 



where: 



The model 2x2 Hamiltonian matrix which was discussed 
above can describe in the NH case two coupled resonance 
states, where Hu is the complex energy of the atomic, molec- 
ular or mesoscopic resonance state that absorbed one photon, 
H\\ = Ei — iTi + tvjj, while #22 = E 2 — ^2 is the complex 
energy of the excited resonance state. T\ > and T2 > are 
the decay rates of the two resonance states. The coupling term 
as usual is given by H\2 = H21 = £0^12/2 where eo is the 
maximum laser field amplitude and di 2 is the complex dipole 
transition matrix element. The NH Hamiltonian matrix can be 
now rewritten such that 



H = 



E 1 + hw+^f 



t<>»12 

2 



2 2 

rp iAT 
Hj2 5- 



ri + r 2 
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1 



(4) 

where AT = Y 2 — T\. As one can see from Eq|4] relative 
gain and loss states are obtained (e.g., when AT > then one 
state has a relative gain while the other resonance state has a 
relative loss). EPis obtained when Eqs. Ill2l are satisfied. Con- 
sequently, an EP in the spectrum is obtained when the maxi- 
mum field amplitude is e^ p = AT / Re[di2] and the laser fre- 
quency is equal to uj ep = {E 2 — Ei — Im[di2]e EP )/h. When 
the laser field is strong enough to allow a multi-photon absorp- 
tion the calculations of the conditions for EP are slightly more 
complicated but achievable. 

In closed systems adiabatic solutions converge to the exact 
solutions of the TDSE in the limit of infinitely slow variation 
of the potential parameters. In open systems the situation is 
very different. In contrast to hermitian systems, in open sys- 
tems for almost any path in parameter space the non-adiabatic 
couplings become more significant the slower the potential 
parameters are varied. As a result the adiabatic theorem often 
breaks down in open quantum systems. Let us explain this for 
our 2x2 model hamiltonian when the two potential param- 
eters qi and q2 are time dependent parameters. These can be 
external field parameters such as laser frequency and intensity. 
The conventional adiabatic approximation is associated with 
the eigenvalues E±(qx, (72) and eigenfunctions 4>'±(qi,q2) of 
the Hamiltonian matrix in Eq|4] The dynamical non-adiabatic 
corrections to the solutions of the TDSE result from the de- 
pendence of the potential parameters qi and q2 on time. The 
dynamical non-adiabatic hamiltonian matrix elements couple 
different adiabatic states and are given by the following matrix 
elements, 



yNA = y +/ _ e +^ T AE ad ( qi (t),q 2 {t))dt 



V + /- = (#. d (9i,«,)|gi£ + &|- \4> ad (qi,q 2 )) 
oq 1 oq 2 

V_ /+ = fa-(9i.Sa)l9io- (6) 
dq 1 dq 2 

and, 

AE ad ( qi , q 2 ) = Ef ( Ql , q 2 ) - El d ( qi , q 2 ) 

= AS ad (q 1 ,q 2 )-iAT ad (q 1 ,q 2 ). (7) 

Here T is the duration of the loop in the parameter space. 
For closed systems where the Hamiltonian is hermitian 
AE ad (gi , q2 ) has real values only. Therefore the exponents 
in Eq[5]are just a phase factor, and the non-adiabatic correc- 
tions vanish as T — > 00 (i.e., the variation of the potential 
parameters is arbitrarily slow). 

In open system where we are dealing with resonances the 
energy difference might be complex. The imaginary part 
of the energy difference, Ar ad > for instance, leads to 
V^ A -> 00 while V^A ->■ as T ->■ 00. The exponential 
divergence in T of the exponent in VffA easily overcomes 



ySA 



V_, + e -ifo ^E ad ( qi (t),q 2 (t))dt 



(5) 
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the 1/T suppression (responsible for the Hermitian adiabatic 
theorem) that is associated with the pre-exponential terms in 
Eq|5]that contain the time derivatives of the potential parame- 
ters qi t 2- This implies that in the limit of slow evolution only 
one state evolves adiabatically while the other state behaves 
non-adiabatically. The adiabatic state is the one which decays 
slower. For other states the adiabatic solution is not valid even 
approximately, making the adiabatic flip often discussed in the 
literature 015U16I1 impossible. 

Let us assume that our protocol implies that the external po- 
tential parameters (qi, 52) are varied in time in a closed loop 
which encircles the EP (see for example Fig|TJ. This EP is 
obtained at the values of the external parameters {qf p 1 5-f )■ 
Note that this dynamical protocol requires to solve the time- 
dependent Scrodinger equation. The only adiabatic solution 
which describes correctly the dynamics is the longest-lived 
resonance state. The lifetimes of the adiabatic resonance 
states are obtained by averaging the decay rates over the en- 
tire closed loop in the potential parameter space. That is, the 
inverse lifetime of the adiabatic states denoted by ± are given 
by T _1 V±dt where is the decay rate at any given 
point on the closed loop. 

A key point in the understanding the difference between the 
dynamics for bound and decaying systems is to realize how 
the non-adiabatic dynamical correction terms couple different 
adiabatic solutions. There is an asymmetric phenomenon in 
the calculation of the strength of the dynamical non-adiabatic 
coupling between two adiabatic resonance solutions which 
do not exist between two bound adiabatic solutions. The 
strength of the non-adiabatic coupling term that induces a 
transition from the \(jf^) state to the \ 4> ad ) state is proportional 
to e~ $0 Ar " ddt while the strength of the coupling that induces 
the transition from the \<fi a ! i ) to the \4>°f) adiabatic resonance 
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FIG. 1 : Two loops in the parameter space of the laser 
intensity and wavelength. When the laser intensity is turned 
on the 9th and the 10th vibrational states of the molecular 
ion become resonances. The EP is obtained when the two 
resonances become degenerate not-hermitian eigenstates. 
Both loops encircle the EP, where the two states coalesce. 



solution is the inverse of this expression, i.e., e + -To Ar dt 
see Eq|5H7] This asymmetric dynamical non-adiabatic effect 
stands behind our discovery that at the end of the propaga- 
tion the system is found in the longest-lived pure state, ir- 
respective of the initial condition. This effect is not due to 
the "evaporation" of the shorter-lived adiabatic states. It hap- 
pens because the short-lived adiabatic states are transformed 
into the long-lived adiabatic states during the time propaga- 
tion process. This phenomenon is first described in our paper 
published in Ref. 11711 . In Ref . 1 1 811 this exchange-state phe- 
nomenon was illustrated for a real physical system when 
interacts with chirped laser pulses. 

We are now at core of the universal asymmetric state ex- 
change phenomenon which is the focus of this study. It is 
based on the observation that the integral J Q T AT ad dt changes 
sign when the loop in parameter space changes from the 
clockwise direction to the counter clockwise direction if and 
only if the closed loop encircles an EP (this is due to the ex- 
change of the instantaneous states which is the property of the 
EP described above). This means that the dynamical proto- 
col described above imposes specific asymmetry for our time 
dependent hamiltonian. The consequences are dramatic. If 
by applying the dynamical protocol in the clockwise direction 
the |</>+ d ) state is obtained as the external parameters return 
to their initial values (independently of the initial state). The 
\4>°^) state will be obtained when the same protocol is applied 
in the counter clockwise direction. This way we can control 
the dynamics and produce a quantum diode-like device for 
atomic, molecular or optical systems, such that the output de- 
pends on the direction in which one enters the device. A gen- 
eral scheme for such device is depicted in Figf2] 

To illustrate and confirm the possibility of a diode-like 
quantum gate by applying the non-adiabatic time-asymmetric 
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(NA-TAS) mechanism we apply the above protocol to 
in a laser field. Fig. [3] demonstrates the interaction of H J 



a|0>+b|1> 



FIG. 2: A schematic representation of a diode-like device 
based on the NA-TAS mechanism. The device can be an 
optical WG (passive diode-like) or an electronic one related 

to any of the applications described below. 
The multi-color line shows the external field as a function of 
the propagation in the device. This field is characterized by 
two parameters, the amplitude and frequency which both 
change during propagation. Two options are presented. In red 
the parameters are varied in a closed loop encircling an EP 
(EPi) whereas in green the EP is outside of the loop (EP2). 
The input in the two cases is the same, i.e. superposition of 
two modes. The output for the two cases is very different. In 
red we show that only when the EP is inside the loop of the 
time varying parameters only ONE mode is obtained, |0) or 
|1). The output depends on the directional propagation in the 
device showing its asymmetric behavior. 



with a chirped laser pulse according to the two different pro- 
tocols shown in Fig. Q] (two different closed loops in laser 
wavelength-intensity parameter plane which encircle the same 
EP). We chose this system since it was a subject of experimen- 
tal studies for many years (see for example Ref. lfl9lo and since 
our theoretical predictions can be experimentally confirmed. 
The numerical results of Fig|3]were obtained by propagating 
the system in the basis of the instantaneous solutions of the 
H 2 interacting with the laser. 

Note that in order to design an experiment we need to add 
the rotational motion into our calculations since the EPs due 
to the coalescence of different rotational states, rather than 
different vibrational states, are more likely to happen. The 
study of the effect of the molecular rotations on the control- 
lable optimal pathway between different ro-vibrational states 
is currently ongoing an is supported by the results presented 
in Fig. [3] 

The results of Fig.|3]can be summarized in Table 1 to illus- 
trate how by going along a loop in one direction we always ob- 
tain one of the coalescing states whereas by going in the other 
direction we will always wind up in the other state. This re- 
sult is independent of the initial superposition of the two field- 
free states. These results clearly demonstrate the power of the 
mechanism proposed in this research with the purpose to in- 
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FIG. 3: The asymmetric state-exchange mechanism in H^. 
W v (t) gives the projection of the propagated wavepacket 
(WP) during a chirped laser pulse on the field-free vibrational 

states of the molecular ions, \v = 9} and \v = 10). In this 
calculation we normalized the propagated WP to unity at any 
given time in order to compensate for the decay. The thin 
black lines describe the solution which starts at the \v = 9) 
vibrational state whereas the thick red lines describe the 
solution which starts at the \v = 10) vibrational state. The 
chirp changes the laser wavelength and intensity in time 
according to the dashed red loop encircling the EP shown in 

Fig. 1. (a) A positive chirped laser pulse corresponds to a 
clockwise trajectory on the loop, (b) A negative chirped laser 
pulse corresponds to a counter-clockwise trajectory on the 
loop. Similar results are obtained when we follow a 
trajectory along the blue solid line in Fig. Q] 
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TABLE I: The driven asymmetric effect on the state-to-state 

vibrational transitions of H2 when either a positive or 
negative chirp laser pulses are used. In a separate column we 
show that the adiabatic approximation yields incorrect 
prediction for the final state. In both cases the laser 
parameters are varied in a closed loop which encircles an EP 
as shown in Fig.Q] 



duce asymmetric transition which is caused by non-adiabatic 
interaction between meta-stable states. 

What is missing in this representation is the fact that the 
NA-TAS mechanism is not only very efficient in producing 
asymmetric molecular diode-like device but also increases the 
quantity of molecules which survive the laser field and have 
not dissociated. In order to illustrate this important property of 
the NA-TAS mechanism we present below the results obtained 
for many different realizations where the duration of the laser 
pulse and the maximum field intensity are varied. We examine 
the effect of the presence of the EP inside or outside the closed 
loop in the chirped laser parameter space on the populations of 
the propagated wavepacket of H^- The parameter p measures 



the proximity of the EP to the loop and when the EP is inside 
the loop then p > while when it is outside p < 0. 
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(a) (b) 

FIG. 4: 

(a) The regions in parameter space where the population of 
the 9th vibrational state of is at least 1000 larger than the 

population of the 10th vibrational state by the end of the 
laser pulse which is positively chirped. Initially at t = the 
two vibrational states of are equally populated. In the 
top panel 10% of the initial population remains, in the middle 
panel 1% remains and in the bottom panel 0.1% survives. 

(b) The regions in parameter space where the population of 
the 10th vibrational state of is at least 1000 larger than 
the population of the 9th vibrational state by the end of the 
lase pulse which is positively chirped. Initially at t = the 
two vibrational states of are equally populated. In the 

top panel 0. 1% of the initial population remains, in the 
middle panel 0.01% remains and in the bottom panel 0.001% 
survives. 



It is evident from the results of Fig J4] that the relevant re- 
gions in parameter space occur for positive p values, i.e., when 
the loop encircles the EP. This clearly shows the efficiency of 
the NA-TAS mechanism first proposed here. The application 
of the NA-TAS protocol results in a diode-like behavior of 
mixed-to-pure state transitions of when the closed loop 
in the frequency-intensity laser parameter space encircles the 
EP. However, above all, the most significant result presented 
here is a fundamental one. It lies in the possibility to observe 
an asymmetric dynamical phenomenon in photo-dissociative 
experiments that reflects the effect of an EP on the dynamics. 
This effect is very different from the effect obtained when the 
degenerate states are bound states and the system is hermitian. 

Here we showed that there are observable physical phe- 
nomena which are hard to explain by hermitian quantum 
mechanics but can be readily explained and predicted by 
using the theoretical tools developed within the framework 
of non-hermitian quantum mechanics. Moreover, the time- 
asymmetric state-exchange mechanism, which is based on our 
ability to locate non-hermitian degeneracies (so called EPs), 
enables one to control the dynamics by external parameters of 
the electromagnetic fields with which the systems under study 
interact. This will open a door to new types of technologies, 
to new type of photo-switches, diode-like atomic, molecular 
and optical devices. 
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